In situ resonant x-ray study of vertical correlation and capping effects during 

GaN/AIN quantum dots growth 



in 
o 
o 

as: 

< 



o 

XJ1 



J. Coraux 1,2 '*, H. Renevier 1,2 , V. Favre-Nicolin 1,2 , G. Renaud 1 , B. Daudin 1 
1 Departement de la Recherche Fondamentale sur la Matiere condensee, 
Service de Physique des Materiaux et Micro structures, 
CEA Grenoble, 17 rue des Martyrs 38054-Grenoble Cedex 9, France 
Universite Joseph Fourier, 
BP 53 F-38041 Grenoble Cedex 9, France 
Johann. Coraux@cea.fr 
(Dated: 2nd February 2008) 

Grazing incidence anomalous x-ray scattering was used to monitor in situ the molecular beam 
epitaxy growth of GaN/AIN quantum dots (QDs). The strain state was studied by means of grazing 
incidence Multi-wavelength Anomalous Diffraction (MAD) in both the QDs and the A1N during the 
progressive coverage of QDs by A1N monolayers. Vertical correlation in the position of the GaN QDs 
was also studied by both grazing incidence MAD and anomalous Grazing Incidence Small Angle 
Scattering (GISAXS) as a function of the number of GaN planes and of the A1N spacer thickness. 
In a regime where the GaN QDs and the A1N capping are mutually strain influenced, a vertical 
correlation in the position of QDs is found with as a side-effect an average increase in the QDs 
width. 
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The control of growth processes is a major issue for the 
achievement of room temperature optoelectronic prop- 
erties in nitride 3D nanostructures such as defect-free 
GaN/AIN quantum dots (QDs) 1 . For that purpose the 
morphological characteristics and structural properties, 
namely the strain, size, size distribution, and density of 
QDs, must be controlled. X-ray techniques have shown 
their powerful abilities to monitor in situ these param- 
eters over large assemblies of QDs 2 . In this letter we 
report on a) an in situ strain analysis of the capping ef- 
fect during the progressive covering of QDs by A1N, car- 
ried out with the Multi-wavelength Anomalous Diffrac- 
tion (MAD) technique in grazing incidence, and b) an 
in situ study of the vertical correlation effect on the 
QDs size during the progressive stacking-up of GaN QDs 
planesiii, carried out by both Grazing Incidence Small 
Angle Scattering (GISAXS) and MAD. 

Samples were grown on 6H-SiC(0001) substrates in 
a plasma assisted molecular beam epitaxy chamber in- 
stalled on the French Collaborating Research Group 
beamline BM32 at the European Synchrotron Radiation 
Facility (Grenoble, France). The nitrogen flux was sup- 
plied by a radio-frequency plasma cell ; Knudsen cells 
provided the Ga and Al fluxes. The ultra-high vac- 
uum environment made possible the in situ growth mon- 
itoring by Reflection High Energy Electron Diffraction 
(RHEED) and x-ray probes. Prior to GaN, a thin 8 
monolayers (MLs) A1N buffer was deposited. The A1N 
surface quality was controlled by RHEED. Wurtzite GaN 
QDs were synthetized on top of A1N via the modified 
Stranski-Krastanow mode&, taking advantage of the in- 
plane 2.5% lattice mismatch between A1N and GaN. The 
QDs were grown at ~ 730° C by depositing a GaN equiv- 
alent amount of 6 MLs. The first sample consisted of 
one GaN QDs plane progressively covered by 2, 4, 6, 8, 
10, 12, 20, 27, 34 MLs of A1N. The second sample was a 
progressive stacking of 10 GaN QDs planes separated by 



25 MLs A1N spacers. The third sample was made of 2 
successive progressive stacks of 5 QDs planes separated 
by a) 25 MLs and b) 53 MLs A1N spacers. 

A grazing incidence and exit setup was used to enhance 
the scattered x-ray signals from the QDs with respect to 
that of the substrate. Diffraction anomalous measure- 
ments were carried out around the Ga K-edge (10.367 
kcV). For the MAD experiments reciprocal space scans 
along the [1010] direction around h=3 were systemati- 
cally recorded at 12 energies around the Ga K-edge us- 
ing a point detector. The detector slits were opened so as 
to measure the integrated intensity over the grazing exit 
angle a/. A 2D charge-coupled device detector placed 
perpendicular to the incident beam was used for GISAXS 
measurements 2 * 

Figure ^a) shows the evolution of /i-scans close to the 
(3030) reflection as a function of the A1N coverage. These 
ft-scans are related to both the in-plane strain state and 
size. The incident x-ray angle was fixed at 0.17°, under 
the critical angle for total external reflection (~ 0.23° 
for bulk GaN in this energy range), so as to limit the 
scattering length below 10 nm and have surface sensitive 
measurements. With no A1N coverage, one observes a 
diffuse-scattering peak ascribed to QDs slightly strained 
by the A1N buffer and SiC substrate. As the A1N cov- 
erage increases, this peak is progressively distorted and 
shifted towards bulk A1N. 

Further analysis was made possible by distinguishing 
the GaN contribution from the A1N and SiC contribu- 
tions, using MAD measurements^^. Figure ^b) shows 
some of the diffracted intensities measured in situ for a 8 
MLs A1N coverage, at 12 energies across the Ga K-edge, 
taking advantage of the Ga anomalous effect to localize 
the Ga contribution along [1010]. For these measure- 
ments, the Ga partial structure factor Fgo, of phase y>Ga, 
that includes the Thomson scattering of all anomalous 
atoms (Ga), can be retrieved. In principle, the retrieval 
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Figure 1: (a) : /i-scans along [lOlO] at 10.267 keV, sub-critical 
incident angle (an = 0.17°), for increasing A1N coverage. The 
SiC substrate peak remains at h=3 while the GaN QDs peak 
progressively shifts to higher h values. Vertical lines indicate 
the bulk GaN and A1N. (b) : Experimental I exp measured at 
4 energies across the Ga K-edge for a 8 MLs A1N coverage, 
(c), (d) : y/I exp measured at 10.267 keV, \\F Ga \\ and \\F N \\ 
extracted for a 8 MLs (c) and a 12 MLs (d) A1N coverage. 



should be run in the framework of the Distorted Wave 
Born ApproximationiSiii. However, one can discard the 
energy dependence at the Ga K-edge of the reflection co- 
efficients, as a consequence of the small Ga amount (6 
equivalent MLs). Therefore, the recorded intensity cor- 
rected for fluorescence, I exv , is proportional to the total 
square structure factor ||F|| 2 . 



i^ll 2 oc 



^exp (-^0 OC 

{[COS (V? T - ^ Ga ) + Pf Ga f + [Sin ((fiT ~ ^Ga) + Pf G cf}( 

where = \\F Ga \\ / (f Ga \\F T \\). The partial structure 
factor Ft of phase <pt that includes the overall contribu- 
tion of non anomalous atoms and the Thomson scatter- 
ing of all anomalous atoms, Faai as well as ipT ~ <PGa, 
can be extracted for all h values, without any structural 
model by fitting eq. ^to the experimental data with the 
NanoMAD algorithmic. 

As shown in fig. ^c,d), Fc a was extracted and Fjv, 
the complex structure factor of non anomalous atoms, 
was deduced from the extracted Ft and lpt — V?Ga, tak- 
ing into account the wurtzite crystallographic structure. 
For h smaller than 2.985, i.e. far from the SiC sub- 
strate peak, only Al and N contribute to F^. The 
h position of the diffuse Fc a and Fpf peak maximum 
is inversely proportional to the in-plane lattice param- 
eter, since (h- h SlC ) /h SlC = (a - a S ic) /asiC where 
asic — 3.08ll and h SlC ~ 3. 

Figure |3 depicts the evolution of the in-plane lattice 
parameters and mismatches in both A1N and GaN, as a 
function of the A1N cap layer thickness. The lattice mis- 
match in the material, relative to the corresponding bulk 
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Figure 2: In-plane lattice parameter and mismatches (relative 
to bulk A1N or GaN) for A1N and GaN deduced from F N and 
FGa extractions. The error bars in the last 3 points for GaN 
are large due to weak extracted Fa a . 
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Figure 3: 2n/ AQ r as a function of the number of QDs planes. 



material, is / = (a — abuik) /dbuik- It is clearly demon- 
strated that the A1N capping initially stressed by the 
QDs (Jain — +1-7%) is then progressively relaxed. Two 
different regimes can be distinguished, in which QDs and 
A1N capping are mutually influenced. A rapid decrease is 
observed until 16 MLs, followed by a more gentle diminu- 
1) tion probably lasting after 34 MLs. Even for large A1N 
deposits (100 nm, not shown in Fig. [5J A1N remains 
slightly expanded (Jain — +0.6%). Jointly the partially 
relaxed (JoaN — —0.7%) surface QDs are gradually com- 
pressed by the A1N capping untill 16 MLs, leading to a 
plateau at ~ —1.1%, keeping in mind the large uncer- 
tainties for A1N coverage thicker than 16 MLs. Interest- 
ingly, the A1N strain state varies only very slowly above 
an A1N thickness of - 30 MLs (7.5 nm). Under this 
semi-quantitative limit, the A1N tickness is such that the 
presence of buried dots significantly affects the A1N strain 
state, and a GaN QDs vertical correlation is expected. 

Another grazing incidence MAD study has been car- 
ried out for a gradual QDs planes stacking. The A1N 
spacer was chosen in a regime of mutual QDs and A1N 
influence and fixed at 25 MLs. Fa a and Fn were ex- 
tracted for 1, 2, 3, 5, 7, and 10 capped QDs planes, in 
order to obtain an information regarding the finite-size of 
GaN domains. Indeed, the full width at half maximum 
AQ r of Fca in radial scans is related to the size Dc a N 
by the Scherrer-like formula 
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Figure 4: (a,b) : GISAXS maps measured at 10.267 keV and 
«j = 0.23° (larger than the critical angle in A1N) for 5 QDs 
planes stackings with (a) 25 MLs and (b) 53 MLs A1N spacers. 
Q z is along the growth direction [0001]. Q x stands for the in- 
plane direction. (c,d) : Experimental (circles) and smoothed 
(lines) GISAXS intensities for Q x = -0.033A -1 (white line 
in (a,c)). (d) is measured at 10.267 keV. The Q z period ob- 
served in (a,c) correspond to 6.8 nm, i.e. ~ the period of the 
stack along (0001). The faint periodic behaviour present in 
(b,d), which period is the same as in (a,c), originates from 
the stacking with 25 MLs A1N spacers underneath. 



where the broadening along the scattering vector 
AQ r = Ah/ (VSasic), Ah being the width of the 
Fca peak along the h direction. Equation assumes 
that the broadening by a strain distribution effect is 
negligible 13 . The evolution of Dc a N with the number 
of plane in the stack is plotted in Figure |31 Despite the 
10% uncertainties, an increase of ~50% in the size of 
QDs is put in evidence, as a consequence of a vertical 
correlation effect. 

The vertical correlation effects were further investi- 



gated by anomalous GISAXS. The incident angle was 
set slightly above the critical angle, making the tech- 
nique sensitive to the whole stacking. Figure E^a) shows 
GISAXS measurements for a stack of 5 QDs planes sep- 
arated by 25 A1N MLs. The occurrence of Q z satellites 
at Q x = ±0.033A _1 is an evidence for the QDs position 
vertical correlation, i.e. along the [0001] directioni^i^. 
Note that the period of these satellites, AQ Z , corresponds 
to a stack period = 2ir/AQ z ~ 6.8nm, i.e. the bilayer 
thickness. The strong energy dependence at the Ga K- 
edge of these satellites (Fig. Etc)) confirms that the scat- 
tering yield is mainly from the QDs. Finally, fig. Bib) 
shows a GISAXS map measured for a QDs planes stack- 
ing with 53 MLs A1N spacers. Compared to the case of 
thiner (25 MLs) A1N spacers (Fig. E(a,c)) 53 MLs A1N 
spacers (Fig. EJb,d)) seem large enough to inhibit the 
vertical correlation effects, likely because they are out 
of the regime of QDs and A1N mutual influence. This 
is consistant with the above conclusion that the critical 
thickness of the A1N spacer requested to allow correlation 
lays around 30 MLs. 

In summary, the in situ diffraction analysis of the cap- 
ping has shown that the strain state for the A1N spacer 
is influenced by the buried QDs up to 30 MLs of A1N 
capping, which corresponds to the upper limit (7.5 nm 
spacers) to observe correlation between successive GaN 
layers. More specifically, the connection between cap- 
ping effects and vertical correlation effects was put in 
evidence. In particular it was shown that in a regime 
where the GaN QDs and the A1N capping are mutually 
strain influenced, a vertical correlation in the position of 
QDs is enabled with, as a side-effect, an average increase 
in the QDs width. 
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